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ABSTRACT 

Context. To better understand the initial conditions of the high-mass star formation process, it is crucial to study at high-angular 
' resolution the morphology, the kinematics, and eventually the interactions of the coldest condensations associated with intermediate- 

I /high-mass star forming regions. 

. Aims. The paper studies the cold condensations in the intermediate-/high-mass proto-cluster IRAS 05345+3157, focusing the attention 

on the interaction with the other objects in the cluster. 

Methods. We have performed millimeter high-angular resolution observations, both in the continuum and several molecular lines, 
with the PdBI and the SMA. In a recent paper, we have already published part of these data. The main finding of that work was the 
^ , detection of two cold and dense gaseous condensations, called N and S (masses ~ 2 and ~ 9 M©), characterised by high values of the 

$_H ' deuterium fractionation (~ 0.1 in both cores) obtained from the column density ratio A'(N2D^)/A'(N2H^). In this paper, we present a 

' full report of the observations, and a complete analysis of the data obtained. 

\ Results. The millimeter maps reveal the presence of 3 cores inside the interferometers primary beam, called CI -a, Cl-b and C2. None 

of them are associated with cores N and S. Cl-b is very likely associated with a newly formed early-B ZAMS star embedded inside 
a hot-core, while Cl-a is more likely associated with a class intermediate-mass protostar. The nature of C2 is unclear. Both Cl-a 
and Cl-b are good candidates as driving sources of a powerful '^CO outflow, which strongly interacts with N, as demonstrated by 
^ ' the velocity gradient of the gas along this condensation. The N2H+ linewidths are between ~ 1 and 2 km s"' in the region where the 

continuum cores are located, and smaller (~ 0.5 - 1.5 km s ') towards N and S, indicating that the gas in the deuterated condensations 
'/^ ■ is more quiescent than that associated with the continuum sources. This is consistent with the fact that they are still in the pre-stellar 

phase and hence the star formation process has not yet taken place in there. 

Conclusions. The study of the gas kinematics across the source indicates a tight interaction between the deuterated condensations 
, and the sources embedded in the millimeter cores. For the nature of N and S, we propose two scenarios: they can be low-mass pre- 

■ stellar condensations or 'seeds' of future high-mass star(s). From these data however it is not possible to establish how the turbulence 

' triggered by the neghbouring cluster of protostars can influence the evolution of the condensations. 

O ■ Key words. Stars: formation - Radio lines: ISM - ISM: individual (IRAS 05345+3157) - ISM: molecules 
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^ 1. Introduction et al. 2005) maintain large abundances in the gas phase, and 

their column density ratios reach values of ~ 0.1 or more, much 

■ - The initial conditions of the star formation process are still higher than the cosmic [D/H] elemental abundance (~ 10 ^ 

a matter of debate. Studies have begun to unveil the chemi- Oliveira et al 2003) 
cal and physical properties of starless low-m ass co res on the 

ve rge of forming low-mass sta i's (Kuiper e t al. 119961: Caselli et The characterisation of the earliest stages of the forma- 

al. I2002al l2002bl; Tafalla et al. l2002l l2006h . demonstrating that tion process of high-mass stars is more difficult than for 

in the dense and cold nuclei of these cores, C-bearing molecular low-mass objects, given their shorter evolutionary timescales, 

sp ecies such as CO and CS are strongly depleted (e.g. Caselli et larger distances, and strong interaction with their environments, 

al. l2002b : Tafalla et al. |2002|) . while N-bearing molecular ions Nevertheless, significant progress has been made by various au- 

such as N2D+ and N2H+ (Caselli et al. 2002a, 2002b; Crapsi thors who have performed extensive studies aimed at the iden- 
tification of massive protostellar candidates, i.e. very young 

~ (< 10^ yr), massive (M> 8M0) stellar objects which have not yet 

Send offprint requests to: F. Fontani, e-mail: • • j j- j- /ivj i- • . i inn/c c -ju 

_ V ^ ■ 1. ionised the surrounding medium (Molinari et al. 1996; Sridharan 

Francesco. Fontani@unige.ch , u„„-i „ . ° , u„„, k ^ i i i i i 

i, „ . , . , ■ , , »»j ™ . J Ti et al. 120021; Fontani et al. I2005D . On the other hand, the pre- 

Based on observations carried out with the IRAM Plateau de Bure — r; — , ' T"^ . ... , , , 

Interferometer. IRAM is supported by INSU/CNRS (France), MPG Protostellar stage, namely the phase in which a cold starless core 

(Germany) and IGN (Spain) evolves towards the onset of the formation of a high-mass pro- 

** The Submillimeter Array is a joint project between the Smithsonian ^ostar, has not yet been investigated in great detail, even though 

Astrophysical Observatory and the Academia Sinica Institute of in the last few years an inc reasing eff^ort has be en dev oted to this 

Astronomy and Astrophysics, and is funded by the Smithsonian study (se e e.g. W ang et al. l2006t Zhang et al. 120091; Beuther & 

Institution and the Academia Sinica. Sridharan |2007|) . 
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Many ultracompact ( UC) H n regions are located in clusters 
(see e.g. Thompson et al. '2006*), and it is common to find mas- 
sive young stellar objects in diff'erent evolutionary stages in the 
same star forming region (Kurtz et al. 2000). Therefore, one ex- 
pects to find even earlier phases of massive star formation in the 
vicinity of UC Hii regions and/or other "signposts" of massive 
star formation activity. With this in mind, several authors have 
recently searched for cold and dense spots in the neighbourhood 
of luminous IRAS point sources (Fontani et al. 2006), UC Hn re - 
gions (Pillai et al. 2007i) and methanol masers (Hill et al. l2006l) . 
In particular, Fontani et al. ([2006) have searched for N2D* emis- 
sion in the gas associated with 10 luminous IRAS sources with 
the IRAM-30m telescope, and detected deuterated gas in 7 out 
of the 10 sources observed. 

The object of the present work is one of the sources stud- 
ied by Fontani et al. ilOO& i. which stands out for its very in- 
teresting characteristics: it is a high-mass (~ ISOM©, Fontani et 
al. l2006) dusty clump, undetected in the MSX Sfim band, located 
nearby the luminous IRAS point source 05345-1-3157 (~ 60" to 
the N-E). The region is located at a distance of 1.8 kpc (Zhang 
et al. 2005), and its surface density (~ 1.3 gr cm"^) and mass- 
to-luminosity ratio (~ SLq/Mq) indicate that it is potentially 
a site of massive star formation (see Fig. 1 of Chakrabarti & 
McKee 12005). Hereafter, we will call our target 1. Previous in- 
terferometric studies have reavealed a clumpy structure of the 
molecular gas i n 1 (M olinari et al. 2002). From IRAM-30m data, 
Fontani et al. (I2OO6I) have measured an average CO depletion 
factor (ratio between expected and observed CO abundance) of 
~ 3 and an average deuterium fractionation (the column den- 
sity ratio between a deuterated species and the corresponding 
one containing hydrogen) of ~ 0.01, three orders of magnitude 
higher than the cosmic [D/H] abundance. These findings would 
indicate for the very first time the possible presence of a high- 
mass p re-stellar core (with T ~ 10 K and «h, ~ lO^cm"-' ), 
analogous to those detected in several low-mass star forming re- 
gions. However, the angular resolution was insufficient to deter- 
mine whether we are dealing with a single high-mass core or 
instead with a sample of low-mass ones. 

Therefore, we have recently mapped 1 at high-angular res- 
olution in the N2H+ (1-0) line with the IRAM Plateau de 
Bure Interferometer (PdBI), and in the N2H+ and N2D+ (3- 
2) lines with the Submillimeter Array (SMA), in order to de- 
rive a detailed map of the deuterium fractionation in the source. 
Simultaneously, we have obtained observations in the continuum 
at ~ 96, ~ 225 and ~ 284 GHz with the two interferometers, 
as well as in several lines of other molecules. The preliminary 
results of these observations, focused mainly on the deuterium 
fractionation derived from the N2H^ and N2D^ data, have been 
published in a recent work (Fontani et al. 2008, hereafter paper 
I). For completeness, in Fig. [T] we summarise the main observa- 
tional results of that work, by showing: (i) the 3 mm continuum 
observed with PdBI, which reveals the presence of 4 cores, one 
of which outside the interferometer primary beam, (ii) the dis- 
tribution of the intensity of the N2H^ (1-0) line, which is ex- 
tended and with a complex structure; (iii) the distribution of the 
N2D^ (3-2) line integrated emission, concentrated in two con- 
densations, called N and S. The integrated emission of the N2H^ 
(3-2) line (not shown in Fig. [TJ is compact, and detected to- 
wards the strongest continuum peak only. In paper I we have 
derived the masses of N and S, which are ~ 9 and ~ 2.5 M0, 
respectively. Also, from the N2D^ /N2H^ column density ratio 
we have obtained a deuterium fractionation of ~ 0.1 in both 
condensations, which are the typical values derived in low-mass 
pre-stellar cores. 



In this work we present a full report of the observations pre- 
sented in paper I, and a complete analysis of the data obtained. 
In Sect. |2] we describe the observations and the data reduction, 
while the results are presented in Sect. [3] In Sect.|4]we derive the 
main physical parameters, which are then discussed in Sect. |5] 
The main findings of this work are summarised in Sect. |6l 



2. Observations and data reduction 

2. 1 . Submillimeter Array observations 

Observations of N2D+ (3-2) (at 231321 MHz) and N2H+ (3-2) 
(at 279511.7 MHz) towards 1 were carried out with the SMA 
(Ho et al. i2004.) in the compact configuration on 30 January 
and 21 February 2007, respectively. The correlator was con- 
figured to observe the continuum emission and several other 
molecular lines simultaneously (see Table [T). The phase centre 
was the nominal position of the sub-mm peak detected with the 
JCMT (Fontani et al. 2006), namely R.A.(J2000)=05''37'"52.4' 
and Dec.(J2000)=32°00'06", and the local standard of rest ve- 
locity Vlsr is -18.4 km s"'. For gain calibrations, observations 
of 1 were alternated with the quasars 3C111 and J0530+135. 
Quasar 3C279 and Callisto were used for passband and flux cal- 
ibration, respectively. The SMA data were calibrated with the 
MIR package (Qi 2005), and imaged with MIRIAD (Sault et 
al. 1995). Channel maps were created with natural weighting, 
attaining a resolution of: 3."7x3."0 for the N2D^ (3-2) channel 
map; 3."0x2."8 for the 225 GHz continuum image; 2"7x2."0 for 
the N2H+ (3-2) channel map; l"9xl."2 for the 284 GHz con- 
tinuum image. At the assumed distance of 1.8 kpc, these values 
translate into spatial resolutions of ~ 0.025 and ~ 0.01 pc at 
225 GHz and 284 GHz, respectively (i.e. 5000 and 2000 A.U., 
respectively). The 3cr level in the 225 and 284 GHz continuum 
images is 0.006 and 0.018 Jy beam"', respectively. 

The spectral resolution in velocity for the N2D+ and N2H^ 
(3-2) was ~0.53 and ~0.44 kms respectively. All the lines 
observed are listed in Table[T] in Cols. 1 and 2 we give the tran- 
sitions observed and the line rest frequency, respectively; Col. 3 
gives the spectral resolution; Col. 4 reports on the detection (Y) 
or non-detection (N) of the transition, and the 3cr level in the 
channel maps of the detected transitions is given in Col. 5. 

2.2. Plateau de Bure Interferometer observations 

We observed the N2H+ (1-0) line at 93173.7725 MHz towards 
1 with the PdBI on 1 1 and 20 August 2006, in the D configu- 
ration, and on 3 April 2007 in the C configuration. In Table [T] 
the spectral resolution in velocity is given. We used the same 
phase reference and Vlsr velocities as for the SMA observa- 
tions. The nearby point sources 0507-1-179 and 0552-1-398 were 
used as phase calibrators, while bandpass and flux scale were 
calibrated from observations of 3C345 and MWC349, respec- 
tively. For continuum measurements, we placed two 320 MHz 
correlator units in the band when making the observations in D 
configuration, and six 320 MHz correlator units in C configu- 
ration. The N2H^ lines were excluded in averaging these units 
to produce the final continuum image (at ~ 96095 MHz). The 
synthesised beam size of the N2H^ channel map was 3."2x3."4, 
while that of the continuum was 3" 1 x3."2 , which translates into 
a spatial resolution of ~ 0.03 pc. The 3cr level in the continuum 
image is 4.2 x lO""* Jy beam"'. We stress that the observations 
of the N2H^ (1-0) and N2D^ (3-2) lines have approximately the 
same angular resolution. The data have been reduced with the 
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Fig. 1. Summary of the main findings of paper I: the dashed contours represent the intensity of the N2H^ (1-0) line integrated 
between -21.2 and -14.5 km s corresponding to the main group of the hyperfine components of this line, observed with the PdBI 
(see also Fig. 2 of paper I). Levels range from the 3cr rms, which is ~ 0.02 Jy beam"', to 0.26 Jy beam"', in steps of 3o". The grey 
scale shows the 96 GHz continuum: the levels range from the 3cr rms (4.2x1 0"'* Jy beam"') to 4x10"-' Jy beam"', in steps of 3cr. 
The solid contours represent the N2D^ (3-2) line emission integrated between -18.37 and -16.7 km s"', obtained with the SMA 
(see also Fig. 1 of paper I). The two main condensations are indicated as N and S. Contour levels start from the 3cr rms (~ 0.09 Jy 
beam"'), and are in steps of 2 cr. The 96 GHz continuum cores are indicated as CI, C2, C3 and C4. The dotted circle represents 
the PdBI primary beam at 96 GHz (~ 48")- The ellipse in the bottom left corner shows the synthesised beam of the N2D^ image, 
comparable to that of the PdBI at 96 GHz. 



GILDAS software, developed at IRAM and the Observatoire de 
Grenoble. 

2.3. Combining SMA and NRAO data for '^CO (2-1 ) 

The problem of the missing short spacing information on the 
SMA data of the '^CO (2-1) line can be solved combining 
the SMA data with the peviously published ' ^CO ( 2-1) map 
from the NRAO- 12m telescope (Zhang et al. l2005h . For ob- 
serving details on the s ingle- dish observations we refer to the 
paper by Zhang et al. ( l2005h . The single-dish data were con- 
verted to visibilities in the MIRIAD package with the task 
UVMODEL. The interferometer and single-dish data were then 
processed together The synthesised beam of the combined data 
is 2.87"x3.60" (PA ~ 42°). For m ore details on the method 
adopted see e.g. Zhang et al. (l2007h . 

2.4. Data reduction 

The N2H+ (1-0), (3-2), and N2D+ (3-2) rotational transitions 
have hyperfine structures. To take this into account, we fitted the 
lines using METHOD HFS of the CLASS program, which is part 
of the GAG software developed at IRAM and the Observatoire 
de Grenoble. This method assumes that all the hyperfine compo- 
nents have the same excitation temperature and width, and that 
their separation is fixed to the laboratory value. The method also 
provides an estimate of the total optical depth of the lines, based 
on the intensity ratio of the different hyperfine components. For 
the rest frequencies of the N2H^ and N2D^ lines, we have used 



the laboratory values given in Crapsi et al. (2005'). We have not 
applied the corrections given in Pagani et al. (2008), since these 
are well below the spectral resolution of our observations. 

3. Results 

3.1. Continuum maps 

The 284 and 225 GHz continuum images obtained with the SMA 
are shown in top and bottom left panels of Fig. |2] respectively. 
For a complete comparison among the three continuum maps, in 
Fig.|2]we also show the continuum emission at 96 GHz, already 
presented in Fig.[T] 

At 284 GHz, three main compact cores are detected inside 
the SMA primary beam (44"), two of them close to the map 
centre and the other one is ~ 10" S-E of the map centre. We 
will call the central condensations Cl-a and Cl-b, respectively, 
because they are located at the position of source CI in the 96 
GHz image, while the eastern one exactly matches the emission 
peak of core C2. Among the other 96 GHz continuum sources 
shown in Fig[T] we clearly detect C3 outside the interferometer 
primary beam (not shown in Fig.|2]l, while C4 is undetected. 

At 225 GHz, the continuum sources Cl-a and Cl-b are 
blended into CI as in the 96 GHz image (Fig. [U, and C2 is 
clearly detected as well. As for the 284 GHz image, C4 is un- 
detected while C3 is detected but not shown because outside the 
primary beam, corresponding to the interferometer field of view. 
There is a faint emission ~ 18" to the NE of Cl-a and Cl-b, 
but we decided not to consider this as a real source because it 
is detected at 6cr in this image only, and quite close to the edge 
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Fig. 2. Top left panel: map of the 284 GHz continuum (grey scale) obtained with the Submillimeter Array. The first level is the 3cr 
rms - 0.018 Jy beam ' (step = 3cr rms). The crosses indicate the position of the 96 GHz continuum sources shown in Fig.[T] except 
C3 which is outside the interferometers primary beam. The thick contours correspond to the integrated emission of the N2D+ (3-2) 
line, starting from the 3cr level and separated by 2cr, and thus indicate the position of the pre-stellar cores detected in paper I (see 
Fig. [TJ. The ellipse in the bottom left corner represents the synthesized beam of the continuum image. Bottom left panel: same as top 
left for the 225 GHz image. The contours start from the 3cr rms (= 0.006 Jy beam"'), and are in steps of 3cr. Right panel: same as 
left panel for the 96 GHz continuum, observed with the PdBI. The levels are the same as in Fig. [1] The dotted contour corresponds 
to the 3cr contour level of the N2H+ (1-0) integrated map shown in Fig. 1. 



of the interferometer primary beam. In both images we also plot 
the 3cr level of the two deuterated cores N and S shown in Fig.[T] 

Because C3 and C4 are outside and on the edge of the pri- 
mary beam, respectively, we decided not to discuss further these 
sources in the following. 

Spitzer Post-Basic Calibrated Data (PBCD) at 24 and 70 yum 
were obtained from the Spitzer Archive. The data were taken in 
Photometry/Super Resolution mode as part of program 20635 
(R. Klein, PI). Observations at 24 microns were obtained on 7 
October 2005 (AORKEY 14944768), with an integration time of 
99 s, and at 70 microns on 2 April 2006 (AORKEY 14945024) 
with an integration time of 76 s (using the Fine Scale mode). 
At present no IRAC observations have been made. The left and 
right panels of Fig. |3]show the Spitzer MIPS continuum images 
at 70 yum and 24 yum. In both images, two peaks roughly coin- 
cident with the submillimeter continuum sources CI and C3 are 
detected. At 70 yum, an extended emission is detected towards 
C2, S and the southern portion of N, but we believe that it is just 
due to the low angular resolution of the map (~ 18"), which do 
not allow to disentangle the contibution of the different sources. 



The 24 yum continuum emission shows two peaks spatially coin- 
cident to the 70yum ones, but they are both more centrally peaked 
than at 70 yum, and have smaller extension, probably due to the 
better angular resolution (~ 6" at 24 yum). A faint emission is 
detected towards C2, at the edge of the main emission peak. No 
emission is detected towards N and S at this wavelength. 



3.2. Molecular line emission. 



In paper I we have shown and discussed the distribution of the 
integrated intensity of the N2D^ and NaH^ lines (see also Fig. [Til. 
However, as said in Sect. 12. 11 several transitions of other species 
have been observed and detected simultaneously to the N2D+ 
and NaH^ ones. In this Section, we complete the analysis started 
in paper I and show a full report of the observations obtained. 
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Transition 



Spec. Resol. Detection 





(GHz) 


(km s-') 




( Jy beam ' ) 




SMA 








(2-1) 


219.560 


0.55 


Y 


0.18 


HNCO (10-9) 


219.798 


0.55 


N 




H2C'30(3i.2-2i,i) 


219.909 


0.55 


N 




SO (5,-4s) 


219.949 


0.55 


Y 


0.12 


CH3OH (80.3-71.3) 


220.079 


0.55 


Y 


0.15 


HCOOCH3 (17-16) 


220.167 


0.55 


N 




"CO (2-1) 


220.399 


0.55 


Y 


0.16 


CH3CN (12-11) 


220.539 - 220.747* 


0.55 


Y 


0.15 


HC,N (83-82) 


220.937 


0.55 


N 




CH3OH (8_,.8-7o.7) 


229.759 


0.53 


Y 


0.18 


CH3OH (32,4-4,,4) 


230.027 


0.53 


N 




'^CO (2-1) 


230.538 


1.0 


Y 


0.15 


OCS (19-18) 


231.062 


0.50 


Y 


0.21 


"CS (5-4) 


231.220 


0.50 


Y 


0.21 


CH3OH (102.2-93,2) 


231.286 


0.50 


N 




NoD+ (3-2) 


231.321 


0.52 


Y 


0.21 


CH36H (9i.,-8o.8) 


278.305 


0.44 


N 




N2H+ (3-2) 


279.512 


0.44 


Y 


1.0 


OCS (23-22) 


279.685 


0.44 


N 




DCO+ (4-3) 


288.144 


0.44 


N 




DCN (4-3) 


288.664 


0.44 


N 




C'^'S (6-5) 


289.209 


0.44 


N 




CH3OH (6-5) 


289.307 - 290.939'' 


0.44 


N 




SiS (16-15) 


290.380 


0.44 


N 




N2H+ (1-0) 


PdBI 

93.173 


0.25 


Y 


0.03 



for the detected transitions only 
*: minimum and maximum frequencies corresponding to the components K=7 and K=0, respectively 
minimum and maximum frequencies corresponding to the components (60,6 - 5o,5) and (6+2,4 - 5+2,3)> respectively 
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Fig. 3. Left panel: map of the 284 GHz continuum (solid white contours) obtained with the Submillimeter AiTay, superimposed on 
the Spitzer MIPS image of 1 at 70 yum, shown by the grey scale (first level and step - 10% of the maximum). The levels of the 284 
GHz image are the same as in Fig.|2] The millimeter continuum sources Cl-a, Cl-b and C2 are labelled as in Fig.|2] The dashed 
contour represents the 3cr level of the N2D^ (3-2) line integrated emission (see Fig.|2]l. Right panel: same as left panel for the 24 
fim Spitzer MIPS image. 



3.2.1 . Integrated emission map of tine opticaily tliin 
component of tlie N2H+ (1-0) line 

The integrated intensity of the main group of hyperfine compo- 
nents of the N2H^ (1-0) line has been shown in Fig. [T] As al- 



ready discussed in Sect.[Tl the emission is extended and presents 
several bumps which are not resolved into separated emission 
peaks. This may be the real shape of the emission in this tracer, 
or it can be due to other reasons, such as an insufficient angular 
resolution, or opacity effects. To better understand this, in Fig.|4] 
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Fig. 4. Integrated intensity of the hyperfine component F\F - 
10 —» 11 of the N2H^ (1-0) line (solid contours), superimposed 
on the 96 GHz continuum (grey-scale). The contours start from 
the 3cr (0.018 Jy beam"'), and are in steps of 3cr. The grey- 
scale is the same as in Fig. |2] The white crosses correspond to 
the three peaks detected in the 284 GHz continuum image, and 
the black crosses mark the position of the NaH^ emission peaks. 
The dashed contrours represent the 3cr level of the N2D^ (3-2) 
integrated emission. The elhpse in the bottom-left corner is the 
synthesised beam. 



we show the distribution of the integrated intensity of the hyper- 
fine component F\F - 10— > 11 only, which is well separated 
from the other components and it is expected to be optically thin: 
from the fitting procedure described in Sect. 12.41 we obtain an 
optical depth for this component of ~ 0.4 ^ 0.7. Fig. |4] shows 
that the extension of the integrated intensity of this component 
is similar to that shown in Fig. [1] but in this case we do detect 
five emission peaks. This suggests that, since the angular reso- 
lution is the same, the lack of well separated intensity peaks in 
Fig- His probably due to optical depth effects. 

The main emission peak is located ~ 4" north of C2, while 
the others are located about ~ 5" east of CI, and ~ 5" south of 
C2, respectively, and two fainter peaks (detected at ~ 6cr) are 
seen towards core N. None of the millimeter sources coincide 
with the five N2H^ emission peaks. These results indicate that 
the source can host other dense cores which are not detected in 
the millimeter continuum observations, probably because they 
are too cold. 



3.2.2. Integrated emission map of tine '^CO (2-1 ) line wings 

Single-dish observations of the sourc e in t he '^CO (2-1) line 
have been performed by Zhang et al. (|2005|) . with a ~ 29 " an- 
gular resolution. They detected a massive outflow in the non- 
Gaussian line wings, the orientation of which is approximately 
in the WE direction, but the angular resolution was not sufficient 
to detect neither the outflow center nor the detailed morphology 
of the outflow itself. The observations reported in this work sig- 
nificantly improve the angular resolution of the previous ones, 
allowing us to better determine the outflow shape and identify 
its origin. 

In Fig. |5] we show the map of the integrated intensity in 
the '^CO (2-1) line wings, derived from channel maps obtained 



combining SMA and NRAO data, as described in Sect. 12.31 The 
blue- and red-shifted emissions have been averaged in the veloc- 
ity intervals (^6.4, -26) km s"' and (-9.2, 6.4) km s"', respec- 
tively. The outflow axis is predominantly oriented in the WE di- 
rection, with redshifted gas in the east and blueshifted gas in the 
west. The lobes are clearly separated and have approximately a 
biconical shape. The outflow center is near the position of the 
continuum source CI, and the exciting source can be either Cl- 
a and Cl-b. The source coincident with the most western peak 
of the N2H^ (1-0) line integrated emission can also contribute 
to the observed emission, but less likely than the continuum 
sources, due to the fact that there are no other signs of proto- 
stellar activity associated with it. The blueshifted gas shows a 
fainter secondary peak to the north of the map. From geometri- 
cal considerations, this northern blue-shifted emission might be 
driven by a source within N, rather than CI -a or Cl-b. However, 
as we will further discuss later, the sources eventually embedded 
within N are probably in the pre-stellar phase, i.e. prior to the 
main accretion phase in which the outflow is expected to form, 
so that this solution seems to us very unlikely. 

The outflow length, from end-to-end, is approximately 35" 
(if we do not consider the secondary peaks of the blueshifted 
emission), corresponding to ~ 0.28 pc at a distance of 1.8 kpc. 
The semi-opening angle is between about 30 and 40°, and the 
spatial separation of the lobes suggests that the inclination angle 
with respect to the line of sight is likely t o be ve ry close to the 
plane of the sky (see also Cabrit & Bertout l2000h . 

In Table |2] we give the line parameters for the calculation 
of the outflow properties, namely: the velocity range of the blue- 
and red- wings, AVb and AVi (Cols. 1 and 2, respectively), the in- 
tegrated intensity of '^CO (2-1) in the wings, J TdV (blue) and 
(red) (Cols. 3 and 4; T is the main beam brigthness temperature), 
and the maximum velocity of the wings, Vmaxt and V'max, (Cols. 5 
and 6), defined as the difference between the maximum velocity 
of the blue- and red- wings and the systemic velocity (i.e. -18.4 
km s"'). The values listed in Cols. 3 and 4 are averaged values 
over the blue and the red lobes, respectively. 



3.2.3. Integrated intensity of other lines 

In Fig.|6]we show the distribution of the integrated intensity of 
all the detected transitions listed in Table [1] averaged over the 
channels with signal, but '^CO (2-1) that has been already dis- 
cussed. 

The integrated intensity of the '^CO (2-1) line peaks at the 
position of CI, but it shows an elongated shape towards the 
north-east, suggesting that some of the emission comes from the 
outflow. On the other hand, the C'**0 (2-1) line emission is less 
extended and it seems to arise from a common gaseous envelope 
in which CI -a and Cl-b are embedded. No significant emission 
of this line is detected towards the '^CO outflow lobes. A similar 
distribution is seen in the integrated emission of '■'CS (5-4) and, 
a bit less, in CH3OH (803-713), showing that these transitions 
trace the same gas. 

The integrated intensity of the two high excitation lines OCS 
(19-18) and CH3CN (12-11) clearly shows the presence of a 
compact hot-core centered on source Cl-b. In both tracers, the 
hot core is unresolved, so that its diameter has an upper limit 
of approximately 3", i.e. ~ 0.026 pc at the given distance. The 
hot-core nature of this condensation is confirmed by the gas tem- 
perature of ~ 200 K that we will derive in Sect. 14.11 

For the CH3OH (8_i,8-7o,7) and SO (56 - 45) line, the emis- 
sion looks clumpy, with a main spot corresponding to the posi- 
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Fig. 5. Top panel: red- and blue-shifted integrated emission of 
the '^CO (2-1) line (combined SMA + NRAO data), superim- 
posed on the 96 GHz continuum map observed with the PdBI 
(grey-scale). Both red and blue contours start from 0.25 Jy 
beam"' (~ 3cr rms), and are in steps of 0.25 Jy beam"'. The 
solid black contours represent the 3cr level of the N2D^ (3-2) 
line integrated emission, and indicate the location of the deuter- 
ated cores N and S (see also Fig.[Tli. The crosses correspond to 
the five N2H^ (1-0) emission peaks shown in Fig.H] The dashed 
circle represents the SMA primary beam at the frequency of the 
'^CO (2-1) line. The beam of the combined '^CO map is shown 
in the bottom left corner Bottom panel: same as top panel but the 
grey-scale corresponds to the continuum at 284 GHz observed 
with SMA. The 3cr level of the N2D+ (3-2) emission and the 
interferometer primary beam are not shown here. 



tion of CI, and another prominent spot coincident with the red- 
shifted '^CO emission, and several fainter spots north and west 
of CI which do not coincide with none of the other tracers ob- 
served. The strong emission of SO (56 - 45) in the red lobe of 



the '^CO outflow is consistent with previous studies, in which 
it has been fo und th at SO abundance is enhanced in outflows 
(Bachiller et al. l200lL Viti et al. l2003l) . and it has been clearly de- 
tected in other massive outflows studied at high-angular reoslu- 
tion (see e.g. the case study of AFGL 5142, Zhang et al. l2007h . 



On the other hand, the CH3OH (S-i s-Toj) line is believed 
to be a Class I methanol maser line, and it has been detected at 
low-angular resoluti on in s everal high-mass star forming regions 
(see e.g. Slysh et al.'2002'). Recently, the SMA has imaged sev- 
eral masers in this line towards the source DR21(OH) (Bourke, 
priv. comm.). Class I methanol masers are often located offset 
from the position of hot cores or HII regions, where Class II 
and other masers, such as those of water and OH, are found. 
They are believed to be coUisionally excited, probably where 
a powerful outflow im pacts the surrounding quiescent material 
(see e.g. Sandell et al. 120031) . However, high-angular resolution 
observations have shown that in some cases they can be found 
in close proximity of HII regions and water masers (Kurtz et 
al. 12004 . The map shown in Fig.|6]clearly show that the (8_i_g- 
Toj) line comes from both the hot core and the red lobe of the 
outflow, indicating that it is originated by shocked gas, and the 
spatial extent is similar to that of SO, even though the SO emis- 
sion is more prominent in the outflow than that of CH3OH. To 
better understand the nature of the CH3OH (8_i,8-7o,7) line, in 
Fig.|2]we show the spectra of this line towards two positions at 
an off'set of (-3", -2") and (-1-5", -1.5") from the map centre. 
These two positions correspond to the main emission peaks of 
the averaged map, and roughly coincide with CI and with the 
inner part of the red lobe of the '^CO outflow, respectively (see 
Fig. |6|. We clearly see that the line in the outflow is much nar- 
rower than the one observed towards core CI (~ 5 km s"' and 
~ 1 kms"', respectively). The broad line seen towards CI is 
likely thermal, while the narrow one observed towards the red 
lobe of the '^CO outflow could come from a single maser spot, 
even though the flux density is not very high (~ 0.08 Jy) and 
could also be thermal as well. 



Finally, to better highlight the emission of '^CO and SO in 
the outflow, in Fig. [8] we show the integrated intensity maps of 
the blue- and red-shifted emission of the '^^CO (2-1) and SO 
(56-45) lines. The blue-shifted emission of '^CO seems to trace 
the cavity created by the blue lobe of the '^CO outflow, while the 
red-shifted emission follows the inner part of the '^CO red lobe 
(see upper panel of Fig. [8]). For SO (56-45), we find a similar re- 
sult for the red-shifted emission, while the blue-shifted emission 
of this line does not follow that of '^CO, but it arises mostly from 
core CI and from the '^CO red lobe. This strange feature could 
be explained by the presence of another outflow not detected in 
'^CO, centered roughly in between core N and core CI. A north- 
ern bump in the blue-shifted emission is also detected, located at 
the position of the northern blue-shifted emission seen in Fig.|5] 
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Fig. 6. Integrated maps of all the molecular Unes listed in Table [T] detected towards i with the SMA. For all transitions, the emission 
has been averaged over all the channels with signal. In each panel, the transition is indicated in the top left corner, the grey scale 
represents the 94 GHz continuum, and the position of the mm cores is indicated by the white crosses. For '^CO, C"^0, CH3OH 
(8_i_8-7o,7) and SO, the contours start from the 3cr level of the averaged map, and are in steps of 3cr, while for the other lines are in 
step of Icr. In the top left panel, where the '^CO (2-1) line is shown, the ellipse in the bottom left corner indicates the synthesised 
beam, and the dashed circle corresponds to the interferometer primary beam. The other lines have the same primary beam and 
similar synthesised beams. 
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Fig. 7. Spectra of the CH30H(8_i,8-7o,7) line obtained towards 
the positions of the two main emission peaks seen in the aver- 
aged map of this line (Fig. |6]). The offsets (in arcsec) are indi- 
cated in the top-right corner, and correspond to the peak position 
of CI and the 'base' of the red lobe of the '^CO outflow. 



4. Derivation of the physical parameters 

4.1. Temperature of core C1-b from CHiCN 

Figure |9] shows the CH3CN (12-11) spectrum at the peak po- 
sition of the emission map shown in Fig. |6] CH3CN is a 
symmetric-top molecule whose rotational levels are described 
by two quantum numbers: J, associated with the total angular 
momentum, and K, its projection on the symmetry axis. Such a 
structure entails that for each radiative transition, J+l — > 7, lines 
with K < J can be seen (for a detailed description see Townes & 
Schawlow 1975). In our observations the bandwidth covers up 
to the K -1 component for the (12-1 1) transition. However, we 
detected only lines up to K=4, for which the energy of the up- 
per level is 184 K. In order to compute the line parameters, we 
have performed Gaussian fits to the observed spectrum assum- 
ing that all the K components arise from the same gas. Hence, 
they have the same LSR velocity and line width. We fixed the 
line separation in the spectrum to the laboratory value and de- 
rive a line width by fitting the spectrum of all the K components. 
Then, under the assumption of LTE conditions, we compute the 
CH3CN spectrum using the gas temperature, the CH3CN den- 
sity, the source size as input parameters, and a line width fixed 
to the observed value (Qiu 2009, private communication). This 
approach does not assume optically thin emission in CH3CN, 
as does in the Boltzmann analysis (see e.g. Kuiper et al. 1984; 
Bergin et al. 1994). The best fit to the CH3CN spectrum shown 
in Fig. |9]yields a range of kinetic temperatures from 150 - 250 
K. These values are comparable to typical temperatures found 
in hot molecular cores surrounding a newly formed massive star 
(see e.g. Kurtz et al. i20()0l) . 



4.2. Physical parameters from dust emission 

The physical parameters derived from the millimeter continuum 
images presented in Sect. 13.11 are listed in Table [3] the position 
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Table 3. Peak position, angular and linear diameter, integrated flux density, mass, H2 volume and column density of the millimeter 
condensations CI (resolved into Cl-a and Cl-b in the 284 GHz image) and C2. The masses are computed for /S - 2, and assuming 
T = 30K. The H2 volume and column densities are calculated assuming a spherical source with diameter equal to the deconvolved 
3cr level. 
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of the emission peak, in R.A. (J2000) and Dec. (J2000), of CI 
(resolved into Cl-a and Cl-b at 284 GHz) and C2 are listed in 
Cols. 2 and 3, respectively; the deconvolved angular (6) and lin- 
ear (D) diameters are given in Cols. 4 and 5. Assuming that the 
emission has a Gaussian profile, we have derived the angular di- 
ameter of the sources deconvolving the observed FWHM with a 
Gaussian corresponding to the synthesised beam of the interfer- 
ometer The observed FWHM has been measured as the geomet- 
rical mean of major and minor axes of the contour at half of the 
intensity peak. The linear diameters have been computed using 
a source distance of 1.8 kpc (Zhang et al. l2005h . 

In Cols. 6 of Table[3]we also give the flux density, Fy, derived 
integrating over the 3cr level in each core. Col. 7 shows the gas 
masses obtained from Fy: assuming constant gas-to-dust ratio, 
optically thin and isothermal conditions, the total gas+dust mass 
is given by: 

Fyd^ 

kyBy{T) 

In Eq. O]), d is the distance, ky is the dust opacity coefficient, 
derived according to ky - ky^(v/vQ'f (where vq - 230 GHz and 
kyg = 0.005 cm^ g ~', which implies a gas-to-dust ratio of 100, 
Kramer et al. ll998h . and By(T) is the Planck function calculated 
at the dust temperature T. For this latter, since we do not have 
a direct estimate (except for Cl-b), we have assumed a 'reason- 
able' value of 30 K (see e.g. Molinari et al. 2000). For Cl-b, 
for which we also have a direct temperature measurement from 
CH3CN (see Sect. 14. It . we give two estimates: one assuming 
r = 30 K and the other assuming the value derived from the 
CH3CN spectrum, i.e. T = 200 K. For the dus t gra in emissivity 
index, we assumed P - 2 (see e.g. Hill et al. 12006 ). From the 
3cr level in the continuum images, we have also derived that we 
are sensitive to point-like sources with masses of 0.6, 0.2 and 
O.OlMo at 284, 225 and 96 GHz, respectively (assuming T = 30 
K). We also estimate a column density sensitivity of ~ 7 x 10^^, 
1.9x 10^^ and 1.4x 10^' cm"^ at the three frequencies mentioned 
above. 



In Cols. 8 and 9 of Table [3] we also list the average H2 
volume and column densities. The H2 average volume densities 
have been computed assuming spherical cores, then the column 
densities have been derived by multiplying the volume densities 
for the core diameters. We have used as angular diameters those 
obtained deconvolving the 3cr levels in the continuum images, 
and not those listed in Table [3] because the continuum fluxes 
have been derived integrating the emission within the 3cr con- 
tours. All volume densities are of the order of 10^ cm"^, while 
the column densities are of the order of 10^^ cm"^, which are 
typical values found in cores associated with intermediate-/high- 
mass forming stars, and are higher than the critical value of 1 gr 
cm"^ required to form a massive star at core centre (Krumholz 
et al.|2008), assuming that the gas is totally made of H2. 

The mass of C2, estimated to be 11.8, 4.7 and ~ 2.9 Mq 
at 96, 225 and 284 GHz, indicates that the core probably hosts 
an intermediate-mass object. The mass of CI is 27 and 11.8 
Mq from the 96 and 225 GHz continuum emission, and when 
it is resolved in two components, Cl-a and Cl-b, in the 284 
GHz image, the masses of the components are 3.8 and 5.2 
Mq, respectively. These values are typical of cores associated 
with intermediate-mass protostellar objects (see e.g. Beltran et 
al. 2008). For Cl-b, adopting the temperature of 200 K from 
CH3CN, we even obtain 0.6 Mq. However, as we will discuss 
better in Sect. 15. II the source embedded in core Cl-b is likely an 
early-B newly formed star. Therefore, the gas mass derived from 
the millimeter continuum represents for this source the mass of 
the circumstellar material only, and not that of the central star. 

Additionally, for all sources these mass estimates, especially 
those derived from the 284 GHz continuum, are very uncertain 
because affected by several problems. First, both the dust tem- 
perature and /3 assumed are very uncertain: /3 is expected to be 
between 1 and 2 in high-mass star forming regions (see e.g. Hill 
et al. 2006 and references therein), while the dust temperature 
can v ary co nsiderably from a source to another (see e.g. Molinari 
et al. l2000t Sridharan et al. 2002). For example, assuming the 
gas temperature derived from NH3 observations, that is T = 17K 



10 



Fontani et al.: The intermediate-/high-mass protocluster in 105345 



1 I I 1 I 1 I I I 

96 GHz ;:<5fif +" '^C0C2-J ) wirfgs 




I ^ I I O . 1 1 I 1 1 1 I 

5'37'"54' 52= 
R.A. (J2000) 

Fig. 8. Top panel: red- and blue-shifted integrated emission of 
the '^CO (2-1) line observed with SMA, superimposed on the 
96 GHz continuum map observed with the PdBl (grey-scale). 
For this latter we use the same contours as in Fig. |5] The in- 
tegration ranges in velocity are (-15.88; -8.14) and (-23.06; - 
20.85) km s"' . The red contours start from 0. 16 Jy beam"' (~ 3cr 
rms), and are in steps of 0.16 Jy beam"'. The blue contours start 
from 0.24 Jy beam"' (~ 3cr rms), and are in steps of 0.16 Jy 
beam"' The solid black contours represent the 3cr level of the 
NaD^ (3-2) line integrated emission as in Fig. |5] The dashed 
circle corresponds to the SMA primary beam at the frequency of 
the '^CO (2-1) line. The red- and blue-dashed lines indicate the 
opening angle of the '^CO outflow lobes as in Fig.|5](excluding 
the northern blue bump). The SMA synthesised beam is shown 
in the bottom left corner Bottom panel: same as top panel for the 
SO (56 - 45) line. The integration ranges in velocity are (-23.5; 
-20.2) and (-14;-3) km s"' . For both red- and blue-shifted emis- 
sion, both the first contour and the step is the 3cr of the averaged 
map (0.06 and 0.09 Jy beam"', respectively). The 3cr level of the 
N2D^ (3-2) emission and the interferometer primary beam are 
not shown here. 



(Jijina et al. ll999h . with = 2 we get 7.9 and 10.9 Mq for Cl- 
a and Cl-b, respectively, and 8.8 and 12.2 assuming T - 17K 
and yS = 1.5 (Mathis & Wiffen ll989|). Second, the 284 GHz 
image is more aff'ected than the others by the problem of the 
missing flux, so that the core masses derived from this image 
are representative of the circumstellar material very close to the 
central object only, while those from the 225 and 96 GHz contin- 
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Fig. 9. Spectrum of the CH3 CN ( 1 2- 1 1 ) fine at the peak position, 
obtained with the SMA. The red line represents the best fit to the 
components K - Qto K = 4, performed as explained in Sect. 14. II 



uum are more representative of the whole envelope. To estimate 
the amount of flux filtered out beacuse of missing short spac- 
ing information, we have compared the flux meas ured w ith the 
single-dish SCUBA map at 850 yum (Fontani et al. l2006h to that 
measured in the 284 GHz map: the peak flux of the single-dish 
map is ~ 2.0 Jy beam"'. Assuming a spectral index of 4, the 
peak intensity at 1 mm (i.e. 284 GHz) is 1.04 Jy beam"'. After 
smoothing the interferometric map at the SCUBA angular reso- 
lution of 14 ", the flux density in the area corresponding to the 
SCUBA beam is 0.11 Jy beam"', which means that we are re- 
covering only ~ 11% of the total flux. Therefore, the masses 
derived from the 284 GHz continuum maps have to be consid- 
ered as lower limits for the circumstellar masses. On the other 
hand, at 96 GHz, the expected flux is 0.01 1 Jy beam"' while the 
value measured from the PdBl map is 0.0075 Jy beam"', so that 
at this frequency we are recovering 68% of the total flux. 

4.3. Physical parameters from NxH^ and NxD^ . 
4.3.1. Deuterium fractionation. 

For completeness, in Table |4] we give the N2H^ and NiD^ total 
column densities, A^(N2H^) and A^(N2D^), and the Deuterium 
fractionation, Dtvac, that we have derived in paper 1. A^(N2H^), 
A^(N2D+) and Dfrac= A?(N2D+)/A^(N2H+), derived for N and S, 
are listed in Cols. 2, 3 and 4, respectively, of the same Table. For 
more details on the methods used as well as in the assumptions 
made, see Sect. 3.1 of paper 1. As already discussed, in both 
condensations, Dfrac is 0. 1 1 , which is comparable to the values of 
Dfiac found in low-mass pre-stellar cores by Crapsi et al. (2005), 
following the same method. These values are also comparable to 
those derived by Filial et al. (2007) in infrared-dark clouds from 
deuterated ammonia. However, their observations are related to 
the very cold, pc-scale molecular envelope, and not to compact 
sub-pc scale cores. 

We could not derive the column densities of N and S from 
the N2H^ (3-2) line because this is undetected towards the two 
molecular condensations, as shown in Fig. 1 of paper 1. Such 
a different distribution of the integrated intensity with respect 
to that of the N2H^ (1-0) line is probably due to the extended 
emission being filtered out differently by the two interferome- 
ters. To compute the amount of the missing flux, we have com- 
pared our interferometric spectra with the single-dish spectra ob- 
tained with the lRAM-30m telescope (see Fontani et al. 2006). 
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Table 4. N2H+ and N2D+ total column densities (A^(N2H+) and 
A^(N2D^)), deuterium fractionation (D^^c), linear diameter (L), 
and mass derived from A^(N2H^) (Mn,h+ ) for the N2D^ conden- 
sations N and S. The uncertainties computed following the stan- 
dard propagation of the errors are given between parentheses. 



source 


A'(N2H+) 


A'(N2D+) 


Ofmc 


L 


Mn,h+ 




(xl0"cm-2) 


(xlO'^cm-^) 




(pc) 


(Mo) 


N 


1.9(0.6) 


2.1(0.3) 


0.11(0.04) 


0.09 


8.7 


S 


1.5(0.5) 


1.6(0.3) 


0.11(0.04) 


0.05 


2.5 




-30 



-20 



-10 



—3 



X 




-25 -20 _y15 
V|sr (km s ) 

Fig. 10. Flux density comparison between IRAM-30m spectra 
and interferometric spectra obtained in the N2H+ (1-0) (top 
panel) and N2H^ (3-2) (bottom panel) lines. The interferomet- 
ric flux density of the N2H^ (1-0) line, measured with the PdBI, 
has been multiplied by a factor of 2, while that of the N2H^ (3- 
2) line, measured with the SMA, by a factor of 5. Both PdBI and 
SMA spectra have been obtained integrating the maps in Figs. 1 
and 2 of paper I over the 3cr level. 



We have resampled the single-dish spectra of N2H^ (1-0) and 
(3-2) to the same resolution in velocity as the interferometric 
spectra. The superposition of the spectra is shown in Fig.[TOl In 
the N2H^ (1-0) line, the flux measured by the PdBI is ~ 2 times 
less than that measured with the IRAM-30m telescope, while in 
the N2H^ (3-2) line the flux measured with the SMA is only one 
fifth of that measured with the 30m antenna. This indicates that 
the extended emission is much more resolved in the SMA map 
of the (3-2) line than in the PdBI map of the (1-0) line. 

4.3.2. Velocity field. 

From the molecular lines observed one can derive information 
about the velocity field in the source. In particular, the peak ve- 



Table 5. Physical parameters of the outflow: column densities 
of the blue- and the red- wing (A^b and A^r), masses (Mb, and 
Mout), momentum (Pout), energy (E), dynamical timescale (fdyn), 
mass entrainment rate (Mout), mechanical force (/"out) and me- 
chanical luminosity (Lm)- All parameters have been calculated 
without correction for the inclination angle (Col. 2), assuming 
the mean inclination angle of 57° (Col. 3) and an extreme incli- 
nation angle of 80° (Col. 4). 
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locity, Vlsr, of the lines gives information on ordered motions 
of the gas (rotation, inward or outward motions), while the line 
widths. Ay provide information on the gas turbulent motions. In 
this section, we show the results obtained from the N2H^ (1- 
0) line, for which the spectra are not affected by central dips as 
those of '^CO and its isotopologues. Also, this tracer is fairly 
well detected towards both the continuum cores and the deuter- 
ated ones, so that it allows us to investigate both ordered and 
turbulent motions across the whole source. 

To derive Vlsr and AV, we have adopted the following pro- 
cedure: from the channel maps of the N2H^ (1-0) line, we have 
extracted a spectrum from each pixel inside the 3cr level of the 
integrated emission (pixel size ~ l"x 1"), and then fitted these 
spectra as described in Sect. 12.41 

In Fig [TT] we show the maps of AV (left panel), and Vlsr 
(right panel) of the N2H^ (1-0) line. From the left panel, one can 
see that the line widths are clearly broader than ~ 1 km s" ' in the 
region where the continuum cores CI and C2 and the N2H^ (1- 
0) main peaks are located. On the other hand, the gas is more 
quiescent around this central region, and towards the position of 
N and S (AV between ~ 0.4 and 1 km s '). These results indicate 
that in the continuum cores, as well as in the N2H^ main peaks, 
the star formation process is actively taking place, while in the 
deuterated cores the gas is more quiescent because they are still 
in the pre-stellar phase. Interestingly, the plot in the right panel 
of Fig. [TT] indicates that in the region where the turbulence is 
higher the gas is also red-shifted with respect to the systemic 
velocity, while no clear trend is seen in the other portions of the 
source, in which the gas velocity is quite close to the systemic 
one. This is probably caused by the interaction with the red lobe 
of the '^CO outflow, which is expanding in that portion of the 
cluster. 



In Sect. 15.21 we will discuss in some more details this inter- 
action, and how AV and Vlsr vary in the deuterated cores N and 
S, by using also the N2D+ (3-2) line. 
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Fig. 11. Left panel: Map of the N2H^ (1-0) line width, AV (coloured-scale). The dotted contours range from 0.4 to 2.0 km s in 
steps of 0.2km s ' .The crosses indicate the position of the continuum sources. The dashed lines represent the 3cr level of N and S. 
The PdBl synthesised beam at the frequency of the N2H+ (1-0) line is shown in the bottom left corner. Right panel: same as top 
panel for the peak velocity, Vlsr, of the N2H'^ (1-0) line. The dotted contours range from -16.8 to -19.2km s in steps of 0.3 
km s"'. 



4.4. Outflow parameters from '^CO line wings 

In Table m the characteristics of the flow are given. Col. 1 lists 
the computed parameters which are: the H2 column density of 
the blue- and red-lobe, N\, and A^i ; the mass of the blue- and 
red-lobe, Mb and Mi, and the total mass, Mout; the momentum. 
Pout, the energy, E; the dynamical timescale, fdyn; the mass en- 
trainment rate, Mout; the mechanical force, Pout; the mechanical 
luminosity, L,^. 

The H2 column densities in both the blue and red lobe, A^b 
and A'r, have been derived from the ^^CO (2-1) line integrated 
emission in the wings according to the standard relations (see 
e.g. Rohlfs & Wilson 2004): 



H2 A^J(b,r) 

CO 



gexp(£'j/rex) 



3k 2J+1 



J(b,r) 



1 - eXp(-T) Jb.r 



^MBCiy . 



(2) 



(3) 



In Eq. 0,^ is the inverse relative abundance of CO (as- 
sumed to be 10"*, Frerking et al. |1982|) . A^j, gj and E] are the 
column density, the statistical weight and the energy (in K) of 
the upper level (i.e. J - 2), respectively, Q is the partition func- 
tion and Tex the excitation temperature. For this latter, we have 
assumed the dust temperature, i.e. 30 K, assuming coupling be- 
tween gas and dust. In Eq. (O, k is the Boltzmann constant, yU the 
dipole moment (0. 1 1 Debye for CO), v the line rest frequency, r 



the optical depth and J TusdV is the integrated intensity in the 
line wings. The values listed in Table |5] have been derived as- 
suming optically thin emission, i.e. ~ 1. By comparing 

the '^CO and '^CO emission in the wings, we have found that 
this assumption is correct for our data. 



All the other parameters given in Table [5 have been derived 
according to Eqs. (2) - (8) of Leurini et al. (l2006h . which take 
into account for correction for the flow inclination angle, 0, with 
respect to the line-of-sight. The parameters are given for no cor- 
rection for the inclination angle, for = 57° (the mean inclina- 
tion angle given by Cabrit & Bertout 2000,), and for an extreme 
value of = 80°. Since the outflow orientation looks close to 
the plane of the sky, we will consider the values computed for 
6 = 57° and - 80° as lower and upper limits, respectively. In 
this range, the flow parameters (outflow mass, momentum, en- 
ergy, mechanical force and mechanical luminosity) are consis- 
tent with the results obtained for outflows associated with high- 
mass YSOs observed at high-angular resolution (e.g. Beuther et 
al. ,2004. .2006.) . 



We will better discuss the nature of the exciting source of the 
CO outflow in Sect.O 
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5. Discussion 

5.1. Nature of the continuum sources 

The maps shown in Sect. 13.11 and the parameters derived in 
Sect. 14.21 allow us to discuss the nature of the millimeter con- 
tinuum cores Cl-b, Cl-a and C2. 

Cl-b probably hosts a very young early-B ZAMS star. This 
is indicated by several pieces of evidence: First, Molinari et 
al. |2002) have detected towards i a cm-continuum source whose 
position roughly corresponds to that of Cl-b, and the Lyman 
continuum derived from their observations suggests that the em- 
bedded object is a newly formed B2-B3 ZAMS star. Second, 
the source is embedded in a hot-core, as demonstrated by the 
kinetic temperature derived from CH3CN (12-11) in Sect. 14.11 
and by the presence of other high-excitation lines (see Fig. |6]l, 
and it lies where the main peak of the mid-infrared 24 fim con- 
tinuum emission is detected (see Fig. O. Third, it is close to 
the center of the outflow lobes detected in '^CO (2-1) (Fig.|5]l, 
even though a significant contribution to the observed blue- and 
red-shifted emission can arise from Cl-a and, less likely, from 
another source not seen in the continuum image, located at the 
position of the western emission peak of N2H+ (see Fig. |5]l. As 
discussed in Sect. 14.21 even though the masses derived form the 
millimeter continuum are just 5.2 and 0.6 Mq, these measure- 
ments represent the mass of the circumstellar gas only, and thus 
they do not contrast with the previous findings. 

The source embedded in Cl-a is probably a very young 
intermediate-mass protostar This is indicated by: the core loca- 
tion, at the centre of the '^CO outflow lobes, which suggests that 
the object embedded in the core is one of the main sources driv- 
ing the powerful outflow; the core gas mass listed in Table[3j the 
non- detecti on of the core in the centimeter continuum (Molinari 
et al. l2002h . indicating an evolutionary stage prior to the forma- 
tion of an HII region. Based on these findings, we believe that 
the source embedded in Cl-a is an intermediate-mass protostar 
still in the main accretion phase (i.e. a Class object). 

The nature of the source embedded in C2 is unclear. The 
core mass (see Table [3]l and the non-detection in the 24 //m im- 
age suggest that it can be an intermediate-mass embedded pro- 
tostar. The source is not associated with a detectable outflow, 
which however could be an indication that the embedded source 
is extremely young. Interestingly, one can notice from Fig. |4] 
that C2 is in the middle of a filamentary structure extended in 
the N-S direction which also contains the deuterated cores N 
and S, and four of the five peaks of the N2H^ (1-0) integrated 
emission. This may indicate that this portion of the source host 
the youngest low- and intermediate-mass objects of the proto- 
cluster, and it is consistent with theoretical models that predict 
that filamentary structures are common morphologies of molec- 
ular clouds i n seve ral classes of dynamical formation models 
(e.g. van Loo 120071) . 

5.1 .1 . A proto-Trapezium system? 

From the positions given in Table |3] we derive that the projected 
linear separation among the three continuum sources goes from 
~ 0.014 pc between Cl-a and Cl-b, to ~ 0.07 pc between Cl- 
a and C2, i.e. from ~ 2800 to 14000 A.U., at the given dis- 
tance. Looking at the 284 and 225 GHz images in Fig. |2] one 
could also argue the possible presence of a fourth continuum 
unresolved source to the NW of Cl-a. This is suggested by the 
elongated shape of Cl-a in this direction, while the synthesised 
beam of the SMA at 284 GHz is elongated in NE - SW direction. 



Assuming spherical symmetry, the stellar density in this region 
is thus ~ 2.3 X 10"* proto-stars pc"-', while it is ~ 1.7. x 10"* 
proto-stars pc"^ assuming the three resolved sources only. These 
values are consistent with the stellar density measured in the 
Trapezium cluster in Orion, estimated to be of the order of 10"* 
stars pc^^*, and derived assuming the lower-mass stars in the clus- 
ter as well. This suggests that the stellar density of ~ 1.7 x 10"* 
proto-stars pc"^ is a lower limit, and that perhaps the three young 
intermediate-/high-mass (proto-)stars embedded into Cl-a, Cl- 
b and C2 (plus possibly the unresolved one to the NW of Cl-a) 
represents a Trapezium-like system in the making. 

Similar multiple systems of forming intermediate-/high- 
mass stars have been discovered recently throu gh high -angular 
resol ution o bservations (see e.g. Beuther et al. l2007bl : Rodon 
et al. l2008h . although these systems are denser than the one 
in 1 (more than 10^ stars pc"^). Neverthless, these stellar den- 
sities are much lower than the value of 10^ stars pc"^ needed 
for the formation of high-mass stars through merging of several 
lower-mass stars at the center of rich clusters (see e.g. Bonnell 
et al. 1998), but comparable to the value suflicient to produce 
binary induced mergers (Bonnell & Bate i2005.) . 

5.2. Nature of the deuterated cores 

In paper 1, we discussed the nature of the deuterated cores based 
on the dust continuum emission and the deuterium fractionation 
derived from the column density ratio A^(N2D+)/A^(N2H+). We 
derived the mass of the cores following two approaches: from 
the virial theorem, and assuming an average abundance of N2H^. 
For completeness, in Col. 6 of Table |4] we list the main results 
obtained from the second approach in paper I. These values, to- 
gether with a Dfrac -0.1 (Col. 4 of Table|4]i, suggested that both 
cores are in the pre-stellar phase but at present they are not mas- 
sive. On the other hand, as discussed in Sect. 3.3 of paper 1, the 
average spectra of the N2H+ (1-0) and N2D"^ (3-2) lines in both 
N and S show lines broader than those typically observed in low- 
mass pre-stellar cores. 

We proposed two possible scenarios for the nature of the 
two N2D+ condensations: (i) they are quiescent low-mass star- 
less cores, and therefore they are going to form low-mass stars 
(in particular, N is likely to harbor several unresolved cores); 
(ii) they are condensations still in dynamic evolution, which can 
be changing their mass and shape, undergoing either fragmenta- 
tion or accretion from the parental cloud. In particular, the sec- 
ond scenario is suggested by the predictions of the 'competitive 
accretion models', which predict formation of massive stars in 
a clustered environment through dynamical interaction among 
several 'seeds' of forming (proto-)stars (see e.g. Zinnecker & 
Yorke 2007 for a review). 

To better understand the nature of the two condensations, a 
detailed analysis of the gas kinematics is certainly very helpful. 
For this reason, here we investigate in more detail the gas kine- 
matics in N and S by using as diagnostics the peak velocity and 
widths of the N2H+ (1-0) and N 2D+ (3 -2) lines, following the 
same approach described in Sect. 14.3.21 

5.2.1. Condensation N 

In the two left panels of Fig. [12] we show the maps of Vlsr and 
Ay, respectively, obtained for N from the N2H^ (1-0) line. The 
same plots derived from N2D+ (3-2) are shown in the two right 
panels. One can see in the first place that the velocity field in 
this core is quite complex. The line width of the N2H^ (1-0) line 
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are slighly larger in the southern region of the core, while that 
of the N2D^ (3-2) looks larger in the northern and, less promi- 
nently, in the south-eastern portion of the core. For both species, 
the values range from ~ 0.6 to ~ 1.4 km s"', and are on aver- 
age 3-4 times larger than those found in low-mass pre-stellar 
cores (see e.g. Caselli et al. 2002a for the case of L1544, Crapsi 
et al. 2005). However, for the N2D^ (3-2) line, in the positions 
where AV ~ 0.5 km s"', these values can be influenced by the 
poor spectral resolution of ~ 0.5 km s ' (see Table[T]i. Therefore, 
it is possible that linewidths smaller than 0.5 km s"' are present 
but not observable with the adopted spectral resolution. The peak 
velocity of both N2H+ and N2D^ lines are between ~ -17.7 and 
~ -18.4 km s"', i.e. close to the systemic velocity, but the gas 
in the southern region of the core looks slightly red-shifted with 
respect to the northern one. Also, from Fig.|5]one can notice that 
the red lobe of the '^CO outflow, likely driven by Cl-b, arises 
from a region spatially coincident with the southern portion of 
core N, and possibly interacts with it: the gas in the red-shifted 
lobe of the outflow probably flows behind the core, dragging 
on part of the southern gas of N. In this way, one can explain 
the observed red-shift of the N2H+ and N2D+ lines in the south- 
ern portion of N. This could also be interpreted as being due to 
the presence of unserolved low-mass cores at a different veloc- 
ity, and the N2H^ emission peaks shown in Fig. |4]is consistent 
with this interpretation. In this case, the observed broad lines 
are just due to the superposition of the different cores. However, 
this seems very unlikely since the red-shifted gas in the southern 
portion of the core is spatially coincident with the bulk of the red 
lobe of the '^CO outflow. So that we believe that the observed 
red-shifted emission and broad lines in the southern portion of 
N are both due to the interaction with the ouflow. 

5.2.2. Condensation S 

In Fig. [13] we show the same maps as in Fig.[T2]for core S. The 
peak velocity of both the N2H+ (1-0) and N2D+ (3-2) lines look 
red-shifted towards the north of the core (which appears to be 
more shifted to the north-east in N2D^), suggestive of a rotation 
motion. A disk-like structure has been observed in the low-mass 
pre-stellar core LI 544, and it is expected to b e caus ed by the 
ambipolar diffiision mechanism (Ciolek & Basu'200d; Caselli et 
al. 2002a). However, L1544 looks also flatter than what is seen 
in core S, which shows a spherical shape. This can be due to 
different orientation with respect to the line of sight and to the 
worse angular resolution of our data, since Caselli et al. (2002a) 
obtained observations with linear resolution a factor of ~ 3 better 
than ours. Interestingly, AV of both N2H^ and N2D^ are larger 
towards the north-eastern part of the core, so that the red-shifted 
gas appears to be also more turbulent. As for N, this is sugges- 
tive of an interaction of the core with the red lobe of the '^CO 
outflow, the southern edge of which touches the upper side of 
core S (see the upper panel of Fig.|5]l. 

In summary, the shape and the mass of this core indicate that 
it is a low-mass pre-stellar core, but the high turbulence and the 
red-shifted gas revealed in the northern part of the core (probably 
triggered by the neighbouring '^CO outflow) can influence the 
evolution of the core itself. 



5.2.3. General conclusion on the nature of N and S 

The main finding of the analysis performed above is that both 
cores are characterised by a probable interaction with the red 



lobe of the CO outflow. This can trigger turbulence in the cores 
themselves and can influence their evolution. In theoretical mod- 
els of clustered star formation, turbulence (which in this case is 
generated by already formed protostars) can create density mod- 
ifications across the cloud, originating several dense and cold 
'seeds', which subsequently can accrete backgound gas that was 
initially not associated with the 'accretion domain' o f the s eed 
itself (see e.g. Bonnefl et al. 2004, and McKee & Tan |2003[) . In 
this scenario, the two condensations can become more massive 
and form pre-stellar cores of higher mass. Alternatively, the in- 
teraction with the other cluster members, in particular with the 
powerful '^CO outflow, could also cause the fragmentation of 
the condensations. 

5.3. Comparison with otiier star forming regions 

5.3.1. Intermediate-/high-mass 

The analysis made in Sects. 15.21 and 15.11 indicates that in 1 
there are objects with masses between ~ 2 and ~ 20 Mq, 
and in different evolutionary stages: pre-stellar core candidates, 
intermediate-mass class protostars, and a newly formed early- 
B ZAMS star. In recent years, many star forming regions harbor- 
ing intermediate- and high-mass young stellar objecs have been 
investigated at high-angular resolution, and several of these re- 
gions show a complex structure as seen in 1. However, only in 
very few of them, evidence of the presence of pre-stellar core 
candidates has been found. 

Among the studies performed so far, 1 appears to have the 
strongest similarities with the high-mass protocluster associated 
with IRAS 20343+4129. Patau et al. (2007) have observed this 
cluster with the SMA in the sub-mm continuum and in '^CO (2- 
1). Since the synthesised beam and the source distance are ~ 3 " 
and 1 .4 kpc, respectively, their observations have linear resolu- 
tion comparable to our ones. They have found two intermediate- 
to high-mass young stellar objects, one of which is driving a 
molecular outflow, and the other one is driving an expanding 
cavity. Interestingly, the gas at the edge of the cavity seems to be 
compressed into several dens e core s, detected in the millimeter 
continuum only. Patau et al. (2001) propose that these are po- 
tential starless cores representing a future generation of stars. In 
this scenario, the formation of new stars in the cluster is caused 
by the interaction between the highest mass objects (which are 
expected to form first at cluster center), and the residual circum- 
stellar (circum-clu ster) m aterial as in 1. 

Leurini et al. (i2007h have investigated at high-angular res- 
olution the high-mass star forming region IRAS 05358+3543, 
and found 4 massive millimeter cores. One of these seems to be 
in the pre-stellar phase, being very cold {T < 20 K) and with a 
large reservoir of material (M ~ 19 M©). However, the core does 
not show any interaction with the other cluster members, and the 
authors suggest that it could be also a very embedded massive 
protostar. 

5.3.2. Low-mass 

Several well-studied low-mass star forming regions show fila- 
mentary structures similar to that seen in 1, and clustered star 
formation. 

In the p-Ophiuchi molecular cloud complex, the large major- 
ity of low-mass millimeter cores associated with both pre- and 
proto-stellar objects are fo und to be d istribu ted into filaments 
(see e.g. Nutter et al. l2006l Andre et al. l2007h . In particular, the 
structure of one of these filaments, Oph A, closely resembles that 
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Fig. 12. Left panels: map of the peak velocity (Vlsr) and line width (AV) derived from the NiH^ (1-0) line inside the 3cr level of 
the N2D^ (3-2) emission of core N. The vertical grey-scale on the right side of the two plots indicates the intensity, in km s of 
Vlsr and AV, respectively. In the bottom left corner, the synthesised beam of the channel maps are shown. Right panels: same as 
left panel for the N2D+ (3-2) Une. 
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Fig. 13. Same as Fig.[T2]for core S. 



of core N. Andre et al. (2007^ have used N2H+ (1-0) to perform 
a detailed study of the kinematics in Oph A, and their results 
are different from ours: first, they derive an average line width 
(thermal + non-thermal component) of ~ 0.4 km s ', which is 
yet more than 2 times smaller than that measured in our deuter- 
ated cores, even though in the NaD^ lines we are limited by a 
spectral resolution of 0.5 km s"', and we are at the limit of our 
resolution towards some parts of the deuterated cores . Second, 
they do not find any evidence of interactions between the differ- 
ent cores, so that there is not feedback from protostellar activity. 
Similar morphologies have been obtained observing N2H^ (1- 
0) towards the Perseus molecular cloud by Kirk et al. (l2007h . in 
which there seem to be a minimal contribution of turbulence in 
the observed line widths. 

These results clearly indicate that pre-stellar cores in low- 
mass star forming regions are commonly found in filaments sim- 
ilarly to what we have found in 1, but they are much less turbulent 
and less interacting with the other cluster members than those 
studied in this work. Therefore, as suggested in Sect. 15. 2[ in 
the condensations associated with 1 the nearby intermediate- and 
high-mass young objects likely play an important role in trigger- 
ing turbulence and, eventually, in causing either the growth or 
the dispersion of the gas in these cores. 

6. Summary and conclusions 

We have presented a full report of the observations partially anal- 
ysed by Fontani et al. (i2008 - paper 1) towards the intermediate- 



/high-mass star forming region IRAS 05345-1-3157. The obser- 
vations have been performed with the PdBI and SMA in the 
following molecular transitions: N2H^ (1-0) with PdBI, N2H^ 
(3-2) and N2D+ (3-2), both with the SMA. At the frequencies 
of these lines, we have simultaneously observed the continuum 
emission, as well as other rotational transitions as '^CO (2-1), 
'^CO (2-1), C"*0 (2-1) and few other high excitation lines of 
less abundant species. In paper I, the main finding presented was 
the detection of two molecular condensations (called N and S) 
showing high values of deuterium fractionation (~ 0.1), derived 
from the column density ratio A^(N2D+)/A^(N2H+). 

In this work, the following main results have been obtained: 

- The continuum emission at 96 and 225 GHz reveals two 
main cores, CI and C2, and CI is resolved into two compo- 
nents in the image at 284 GHz, called CI -a and Cl-b. CI is 
coincident with a prominent mid-infrared emission detetcted 
in the Spitzer MIPS images at 24 and 70 /urn, which is only 
barely detected towards C2. 

- The integrated intensity of the optically thin component of 
the N2H^ (1-0) line shows an extended distribution with 5 
emission peaks. Two of them fall inside condensation N, an- 
other one at the edge of condensation S, while none of them 
overlaps with the continuum sources. 

- The integrated emission in the '^CO (2-1) line wings reveals 
the presence of a powerful bipolar outflow oriented roughly 
in the WE direction. From the outflow geometry, the sources 
CI -a and Cl-b are the best candidates for powering the out- 
flow. Assuming a unique source driving the outflow, its pa- 
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rameters are consistent with an engine which is a high-mass 
young stellar object. 

- The N2H^ (1-0) line widths are between 1 and 2 km s ' in 
a region that spatially corresponds to that where the contin- 
uum cores are located, while they are significantly smaller 
(between 0.5 and 1.5 km s"') in the deuterated cores. 

- Based on previous observations and on the results presented 
in this work, we can conclude that Cl-b very likely harbors 
a newly formed early-B star embedded inside a hot-core, 
and CI -a harbors an intermediate-mass class protostar The 
nature of C2 is unclear but it could be a very embedded 
intermediate-mass protostar If we consider the system of 
these three continuum sources and roughly assume a spher- 
ical symmetry, we deduce a star density of ~ lO'* stars pc 

consistent with the stellar density of the Trapezium clus- 
ter and much lower than the value required to form massive 
stars through collisions of low-mass ones. 

- The nature of the deuterated cores has been discussed using 
mainly the information on the gas kinematics that one can 
derive from the lines of N2H^ and N2D^. Core S is likely 
a single low-mass pre-stellar core, and the velocity of both 
the NaH^ and NaD^ lines indicate a rotation motion roughly 
in the N-S direction. This motion and the high turbulence 
measured at the northern edge of the core is likely due to 
dynamical interaction with the red lobe of the '^CO powerful 
outflow. The nature of core N is less clear Two scenarios are 
possible: it can be a filament of unresolved low-mass pre- 
stellar cores, similarly to Oph A, and in this case the high 
values of the N2H^ and N2D^ line widths are just due to the 
superposition of the several unresolved cores. Alternatively, 
it can be a single object still in dynamic evolution, which 
actively interact with the red lobe of the '^CO outflow. In this 
case, the triggered trubulence can cause either the growth or 
the fragmentation of the core itself. 

From these observations, it is not at all clear how the interac- 
tion with the other cluster members can make the condensations 
evolve towards the formation of low- or high-mass objects. Only 
higher resolution kinematic studies of the large scale gas around 
these condensations, and the comparison of these results with 
hydrodynamic evolution models, will help to understand the core 
evolution and provide important constraints on the dynamical 
evolution of intermediate-/high-mass star forming regions. 
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